A chaotic system with only one stable equilibrium 
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If you are given a simple three-dimensional autonomous quadratic system that has only 
one stable equilibrium, what would you predict its dynamics to be, stable or periodic? Will 
it be surprising if you are shown that such a system is actually chaotic? Although chaos 
theory for three-dimensional autonomous systems has been intensively and extensively 
studied since the time of Lorenz in the 1960s, and the theory has become quite mature 
today, it seems that no one would anticipate a possibility of finding a three-dimensional 
autonomous quadratic chaotic system with only one stable equilibrium. The discovery of the 
new system, to be reported in this Letter, is indeed striking because for a three-dimensional 
autonomous quadratic system with a single stable node-focus equilibrium, one typically 
would anticipate non-chaotic and even asymptotically converging behaviors. Although the 
new system is of non-hyperbolic type, therefore the familiar Si'lnikov homoclinic criterion 
is not applicable, it is demonstrated to be chaotic in the sense of having a positive largest 
Lyapunov exponent, a fractional dimension, a continuous broad frequency spectrum, and a 
period-doubling route to chaos. 

PACS: 05.45.-a, 05.45.Ac, 05.45.Pq 
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I. INTRODUCTION 

For three-dimensional (3D) autonomous hyperbolic type of chaotic systems, a commonly ac- 
cepted criterion for proving the existence of chaos is due to Si'lnikov [1-4], which has a slight 
extension recently [5]. Chaos in the Si'lnikov type of 3D autonomous quadratic dynamical systems 
may be classified into four subclasses [6]: 

• chaos of the Si'lnikov homoclinic-orbit type; 

• chaos of the Si'lnikov heteroclinic-orbit type; 

• chaos of the hybrid type with both Si'lnikov homoclinic and heteroclinic orbits; 

• chaos of other types. 
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In this classification, a system is required to have a saddle-focus type of equilibrium, which 
belongs to the hyperbolic type at large. 

Notice that although most chaotic systems are of hyperbolic type, there are still many others 
that are not so. For non-hyperbolic type of chaos, saddle-focus equilibrium typically does not exist 
in the systems, as can be seen from Table I which includes several non-hyperbolic chaotic systems 
found by Sprott [7-10]. More recently, Yang and Chen also found a chaotic system with one saddle 
and two stable node-foci [11] and, moreover, an unusual 3D autonomous quadratic Lorenz-like 
chaotic system with only two stable node-foci [12]. In fact, similar examples can be easily found 
from the literature. 

In this paper, we report a very surprising finding of a simple 3D autonomous chaotic system 
that has only one equilibrium and, furthermore, this equilibrium is a stable node-focus. For such 
a system, one almost surely would expect asymptotically convergent behaviors or, at best, would 
not anticipate chaos per sc. 

From Table I, one may observe that the Sprott D and E systems also have only one equilibrium, 
but nevertheless this equilibrium is not stable. From this point of view, it is easy to understand 
and indeed easy to prove that the new system will not be topologically equivalent to the Sprott 
systems. 

II. THE NEW SYSTEM 
A. The mechanism of generating the new system 

The mechanism of generating the new system is simple and intuitive. 

To start with, let us first review some of the Sprott chaotic systems listed in Table I, namely 
those with only one equilibrium. One can easily see that systems I, J, L, N and R all have only one 
saddle-focus equilibrium, while systems D and E both degenerate in the sense that their Jacobian 
eigenvalues at the equilibria consist of one conjugate pair of pure imaginary numbers and one real 
number. Clearly, the corresponding equilibria are not stable. 

It is also easy to imagine that a tiny perturbation to the system may be able to change such a 
degenerate equilibrium to a stable one. Therefore, we added a simple constant control parameter 
to an aforementioned Sprott chaotic system, trying to change the stability of its single equilibrium 
to a stable one while preserving its chaotic dynamics. 
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TABLE I: Equilibria and eigenvalues of several typical Sprott systems. 
Systems Equations Equilibria Eigenvalues 



Sprott x = -y (0, 0, 0) 0, ±i 

Case D y = x + z 

z = xz + 'iy^ 

Sprott x = yz (0.25,0.0625,0) -l,±0.5i 

Case E y = x"^ —y 
z = l- 4x 

Sprott x = -0.2y (0,0,0) -1.13449, 0.06725 ± 0.58996i 

Case I y = x + z 

z = x + y'^ — z 

Sprott x = 2z (0,0,0) -2.31460, 0.15730 ± 1.30515i 

Case J y = —2y + z 
z = -X + y + y'^ 

Sprott x = y + 3.9z (1,0.9,-0.23077) -1.43329, 0.21664 ± 1.63526i 
Case L y = 0.9x^ - y 
z = 1 — X 

Sprott x = -2y (-0.25,0,0.5) -2.31460, 0.15730 ± 1.30515i 

Case N y = x + z^ 
z = l + y-2z 

Sprott X = 0.9-2/ (-0.44444, 0.9, -0.4) -1.23212, 0.11606 ± 0.84674i 
Case R y = 0.4 + ^; 

z = xy — z 

As a result, we obtained the following new system: 

X = yz + a 

' ij = x^ -y (1) 
z = 1 — 4x. 

V. 
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When a = 0, it is the Sprott E system; when a ^ 0, however, the stabihty of the single 
equihbrium is fundamentally different, as can be verified and compared between the results shown 
in Table I and Table II, respectively. 

TABLE II: Equilibria and eigenvalues of the new system. 
Systems Equations Equilibria Eigenvalues 

New System i = + a (0.25,0.0625,0.08) -1.03140, 0.01570 ± 0.49208i 
a=-0.005 y — — y 

z = 1 - Ax 

New System x = yz + a (0.25, 0.0625, -0.096) -0.96069, -0.01966 ± 0.50975i 
a=:0.006 y = x^ — y 

z = 1 - ix 

New System x = yz + a (0.25, 0.0625, -0.352) -0.84580, -0.07710 ± 0.53818i 
a=0.022 y = x'^-y 

z = 1 - ix 

New System i + a (0.25,0.0625,-0.48) -0.78217, -0.10891 ± 0.55476i 
a=0.030 y = x^-y 

z = 1 — Ax 

New System i = + a (0.25,0.0625,-0.8) -0.60746, -0.19627 ± 0.61076i 
a=0.050 y = x'^-y 
z = 1 — Ax 

To better understand the new system ([T]), and more importantly to demonstrate that this new 
system is indeed chaotic, some basic properties of the system are briefly analyzed next. 

B. Equilibrium and stability 

The system ([1]) possesses only one equilibrium: 

P(xB,y£;,Zij) = Q,^,-16aj . (2) 
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Linearizing the system at the equihbrium P gives the Jacobian matrix 



J 



z y 
2x -1 
-4 





-4 

By solving the characteristic equation |A/ — J| = 0, one obtains the Jacobian eigenvalues, as shown 







-16a ^ 



-1 



(3) 



in Table II for some chosen values of the parameter a. 



C. Lyapunov exponents 

To verify the chaoticity of system ([T|), its Lyapunov exponents and Lyapunov dimension are 
calculated. 

The Lyapunov exponents are denoted by Li, i = 1,2, 3, and ordered as Li > L2 > L3. A system 
is considered chaotic if Li > 0, L2 = 0, L3 < with \Li\ < {L^l. 
The Lyapunov dimension is defined by 

1 ^ 
l^i+il 

where j is the largest integer satisfying X^^^]^ Li > and X]i=i < 0. 

FIG. [1] shows the dependence of the largest Lyapunov exponent of system ([T|) on the parameter 
a. From FIG. [H it is clear that the largest Lyapunov exponent decreases as the parameter a 
increases from —0.01 to 0.05. 



D. The degenerate case of a = (Sprott E system) 

When a = —0.005, the system equilibrium is of the regular saddle-focus type; this case of the 
chaotic system has been studied before therefore will not be discussed here. 

When a = 0, the equilibrium degenerates. It is precisely the Sprott E system listed in Table 
I (see FIG. [2]). The Si'lnikov homoclinic criterion might be applied to this system to show the 
existence of chaos, however, but it involves somewhat subtle mathematical arguments. 

In this degenerate case, the positive largest Lyapunov exponent of the system (see Table II) 
still indicates the existence of chaos. In the time domain, FIG. [3] (top part) shows an apparently 
chaotic waveform of y{t); while in the frequency domain, FIG. [3] (bottom part) shows an apparently 
continuous broadband spectrum |y(i)|. These all prove that the Sprott E system, or the new system 
([1]) with a = 0, is indeed chaotic. 
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FIG. 1: (Color online) The largest Lyapunov exponent versus the parameter a. 



Three-dimensional view x-y pinase plane 




-1 1 2 0.5 1 1.5 2 

X Y 



FIG. 2: (Color online) The new system: chaotic attractor with a = 0, including 3D views on the x-y plane, 
x-z plane and y-z plane. 
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Time series of y(t) 
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FIG. 3: (Color online) Top: An apparently chaotic waveform of y{t) (a = 0). Bottom: An apparently 
continuous broadband frequency spectrum |y(i)|. 

E. The case of a = 0.006: a new type of chaos 



When a > 0, the stability of the equilibrium is fundamentally different from that of the Sprott E 
system. In this case, the equilibrium becomes a node-focus (see Table II). The Si'lnikov homoclinic 
criterion is therefore inapplicable to this case. 

Take a = 0.006 as an example. Numerical calculation of the Lyapunov exponents gives Li = 
0.0489, L2 = and L3 = —1.0485, indicating the existence of chaos. 

In the time domain, FIG. [5] (top part) shows an apparently chaotic waveform y{t); while in 
the frequency domain, FIG. [5] (bottom part) shows an apparently continuous broadband spectrum 
|y(t)|. These all prove that the new system ([1]) with a = 0.006 is indeed chaotic. 
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Three-dimensional view x-y phase plane 
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FIG. 4: (Color online) The new system: chaotic attractor with a — 0.006, including 3D views on the x-y 
plane, x-z plane and y-z plane. 

F. Bifurcations analysis 

FIG. [6] shows a bifurcation diagram versus the parameter a, demonstrating a period-doubling 
route to chaos. 

FIG. [7] also demonstrates the gradual evolving dynamical process as a is continuously varied. 
Both figures indicate that although the equilibrium is changed from an unstable saddle-focus 
to a stable node-focus, the chaotic dynamics survive in a relative narrow range of the parameter a. 
All the above numerical results are summarized in Table III. 

III. CONCLUSION 



This paper has reported the finding of a simple three-dimensional autonomous chaotic system 
which, very surprisingly, has only one stable node-focus equilibrium. The discovery of this new 
system is striking, because with a single stable equilibrium in a 3D autonomous quadratic system. 
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Time series of y(t) 
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FIG. 5: (Color online) Top: An apparently chaotic waveform of y{t) (a = 0.006). Bottom: An apparently 
continuous broadband frequency spectrum |y(i)|. 

one typically would anticipate non-chaotic and even asymptotically converging behaviors. Yet, 
unexpectedly, this system is chaotic. Although the new system is non-hyperbolic type, therefore 
the Si'lnikov homoclinic criterion is not applicable, it has been verified to be chaotic in the sense 
of having a positive largest Lyapunov exponent, a fractional dimension, a continuous frequency 
spectrum, and a period-doubling route to chaos. 

Although the fundamental chaos theory for autonomous dynamical systems seems to have 
reached its maturity today, our finding reveals some new mysterious features of chaos. 
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FIG. 6: (Color online) Bifurcation diagram, showing a period-doubling route to chaos in y (at x = 0.25) 
versus the parameter a. 
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TABLE III: Numerical results for some values of the parameter a with initial values (1,1,1). 
Parameters Eigenvalues Lyapunov Exponents Fractal Dimensions 



a = -0.005 Ai = -1.03140 Li = 0.0884 

A2,3 = 0.01570 ±0.49208z ^2=0 L>l = 2.081 

L3 = -1.0884 

a = Ai = -1 Li= 0.0766 

A2,3 = ±0.5i L2 = Dl = 2.071 
L3 = -1.0766 

a = 0.006 Ai = -0.96069 Li = 0.0510 

A2,3 = -0.01966 ± 0.50975i L2 = Dl = 2.048 

L3 = -1.0510 

a = 0.022 Ai = -0.84580 Li = 

A2,3 = -0.07710 ± 0.53818i L2 = -0.1381 Dl = 1.000 

L3 = -0.8619 

a = 0.030 Ai = -0.78217 Li = 

A2,3 = -0.10891 ± 0.55476i L2 = -0.0826 Dl = 1.000 

L:i = -0.9174 



a = 0.050 Ai = -0.60746 Li = 

A2,3 = -0.19627 ± 0.61076i L2 = -0.0518 Dl = 1.001 
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(c) (d) 
FIG. 7: (Color online) Phase portraits and frequency spectrums: (a) a — 0.006, (b) a — 0.022, (c) a = 0.03, 
(d) a = 0.05. 



